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Abstract 

In this talk, after a short phenomenological introduction on glasses, I will describe 
some recent progresses that have been done in glasses using the replica method in 
the definition and in the evaluation of the configurational entropy (or complexity). 
These results are at the basis of some analytic computations of the thermodynamic 
glass transition and of the properties below the phase transition point. 
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1 Introduction 



Glasses are roughly speaking liquids that do not crystallize also at very low 
temperature (to be precise: glasses also do not quasi-crystallize). 

These liquids can avoid crystalization mainly for two reasons: 

• The liquid does not crystallize because it is cooled very fast: the crystallisa- 
tion time may become very large at low temperature (e.g. hard spheres at 
high pressure). The system should be cooled very fast at temperatures near 
the melting point; however if crystalization is avoided, and the temperature 
is low enough, (e.g. near the glass transition) the system may be cooled very 
slowly without producing crystalization. 

• The liquid does not crystallize even at equilibrium. An example is a binary 
mixture of hard spheres with different radius: 50% type A (radius ta, 50% 
type B (radius r% where R denotes r B /r A . If .77 < R < .89 (the bounds 
may be not precise), the amorphous packing is more dense than a periodic 
packing, distorted by defects. 

Which of the two mechanism is present is irrelevant for understanding the 
liquid glass transition. 
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Fig. 1. The viscosity as function of the temperature according to the Vogel Fulcher 
law (Tjc = 1 full line) and to the mode coupling theory (Tk = 1.375 dotted line). 

Other examples of glassy systems are mixtures: in this case we have M kind 
of particles and the Hamiltonian of N particles is given by 

h= E E E v a>b (x a (i)- Xb (k)) , (i) 

a,b=l,M i=l,N(a) k=l,N(b) 



where a = 1, M; N(a) = Nc(a); J2 a c i a ) = 1- m this case the values of the N 
concentrations c and the MiM + l)/2 functions V a , b (x) describe the model. 
Well studied case are: 

• Soft spheres (non-realistic, but simple), e.g. M — 2, c(l) = .5, c(2) = .5, 
V ayb {x) = R a ,bX~ 12 - 

• Lennard- Jones spheres (more realistic), e.g. M — 2, c(l) = .8, c(2) = .2, 

V a>b (x) = Ra,bX~ 12 ~ A a:b X~ 6 . 

There are some choices of the parameters which are well studied and in these 
case it is known that the system does not crystallize, one of these has been in- 
troduced by Kob and Anderson in the L-J case and correspond to a particular 
choice of the parameters R and A. 

There are many other material which are glass forming, e.g. short polymers, 
asymmetric molecules (e.g. OTP) . . . 

The behaviour of the viscosity in glass forming liquid is very interesting. We 
recall that the viscosity can be defined microscopically as follows. We consider 
a system in a box of large volume V and we call T^ v (t) the total stress tensor 
at time t. We define the correlation function of the stress tensor at different 
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Fig. 2. The specific heat and the entropy excess AS for various cooling rates, 
times: 

(T^(t)T PiCT (0)> = VS mu {t) (2) 

One finds, neglecting indices, that rj oc J dtS(t) ~ r~ Q , where r is the char- 
acteristic time of the system (in a first approximation we can suppose that 
a = 1). 

In fragile glasses the behaviour of the viscosity as function of the temperature 
can be phenomenological described by the two following regimes: 

• In an high temperature region the mode coupling theory is valid: it predicts 
7] oc (T — T c )~ 7 , where 7 is not an universal quantity and it a of the orders 
of a few units. 

• In the low temperature region by the Vogel Fulcher law is satisfied: it pre- 
dicts that rj oc exp(A(T — T^) _1 ). 

Nearly tautologically fragile glasses can be defined as those glasses that have 
Tk 7^ 0; strong glasses have Tk ~ 0. 

The glass temperature (T g ) is the temperature at which the viscosity becomes 
so large that it cannot be any more measured. This happens after an increase of 
about 18 order of magnitude (which correspond to a microscopic time chang- 
ing from 10 -14 to 10 4 seconds): the relaxation time becomes larger than the 
experimental tine. 

A characteristic of glasses is the dependance of the specific heat on the cooling 
rate. There is a (slightly rounded) discontinuity in the specific heat that it is 
shifted at lower temperatures when we increase the cooling time. 

For systems that do crystallize if cooled very slowly, one can plot AS = 
S{liquid) — S(crystal) versus T in the thermalized region. One gets a very 
smooth curve whose extrapolation becomes negative at a finite temperature. 
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A negative AS does not make sense, so there is a wide spread belief that 
a phase transition is present before (and quite likely near) the point where 
the entropy becomes negative. Such a thermodynamic transition (suggested 
by Kaufmann) would be characterized by a jump of the specific heat. Quite 
often this temperature is very similar to the temperature found by fitting the 
viscosity by the Volker Fulcher law and the two temperature are believed to 
coincide. 

Generally speaking it is expected that a thermodynamic phase transition in- 
duces a divergent correlation time. However the exponential dependance of the 
correlation time on the temperature is not so common (in conventional criti- 
cal slowing down we should have a power like behaviour); moreover an other 
strange property of the glass transition is the apparent absence of a detectable 
correlation length or susceptibility (linear or non-linear) which diverges when 
we approach the transition point. 

In section 2, I will firstly present the many valley picture which has been the 
inspiration point of many works; in the subsequent section I will summarize 
some of the results that have been obtained using the replica theory on the 
organization of these valley in phase space, on the temperature dependence of 
the free energy in each valley, in brief of the free energy landscape. A strategy 
for applying these ideas to glasses will be presented in section 4 and in the 
next section I will show some results that have been obtained in the case of a 
Lennard- Jones binary mixture. 



2 The many valleys picture 



A quite old statement is that the glass is frozen liquid, which I prefer to 
reformulate it as that a low temperature a liquid is an unfrozen glass. 

The general idea is the following: the phase space of the systems can be approx- 
imately decomposed into valleys, labelled by a, such that the barrier among 
these valleys is large in the liquid (in the low temperature region where the 
viscosity is large) the barriers become infinite at T K [1]. 

For simplicity I will assume that valleys do maintain their identity when chang- 
ing the temperature, (in reality they split into smaller valleys when we decrease 
the temperature [2]). Under this assumption there is a one to one correspon- 
dence among valleys and inherent structures, i.e. minima of the Hamiltonian 
H . The properties of the liquid are therefore connected to the properties of 
the system near the minima of the Hamiltonian. 
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Fig. 3. An artistic view of a corrugated free energy landscape, the horizontal axis 
being the configuration space. 

The partition function can be approximately written as 



where f a {/3) is the free energy density of the a th valley. As we shall see later 
in the infinite volume limit the previous sum is dominated by those valleys 
having a given free energy density / = /*(/?). 

The number of relevant valleys at / = /*(/?) (i.e. Af*(/3)) is supposed to 
exponentially increase with the size (N) of the system: 



It can be argued that the configurational entropy or complexity £*(/?) is ap- 
proximately near to AS: indeed it has suggest long time ago that the config- 
urational entropy £*, goes to zero at the phase transition point [1]. 

Recently using the techniques stemming from replica theory we have added 
new ingredients to this old scenario: 

• We have construct a microscopic realization of the previous ideas in the 
mean field models [3,4]. In other words there are soluble infinite range mod- 
els in which there is an exponentially large number of valleys and their 
properties can be computed analytically. 

• If we assume that minimal corrections are present in finite dimensions to 




(3) 



a 



cN* ^exp(NE*(/3)) . 



(4) 
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mean field predictions, the mean field results may be extended to three 
dimensional glasses. 
• The replica method [5] gives tools to do the appropriate computations, both 
in mean the mean field approach and in short range models. 

As outcome of these advances a many valleys picture with detailed predictions 
on the landscape has been constructed and the properties of glasses can be 
computed in the framework of replicated liquid theory [6-9]. 



3 The free energy landscape 



Let us suppose that we can define a temperature free energy functional 
as function of the density p(x). It is natural to assume that the different 
valleys are in a first approximation associated to different local minima of this 
functional, the free energy evaluated at the local minima being the free energy 
of the corresponding valley. 

I will summarize which are in mean field the properties of the local minima the 
free energy that have been computed explicitly in some long range microscopic 
models. Although these results have been obtained in some particular models, 
one can show that they are valid for a quite large class of models (at least in 
the mean field approximation [10,11]). 

• For T > T c there is one relevant minimum with p(x) = const. There are 
many solutions of the equation 

which are not minima (there are some negative eigenvalues of the Hessian 

• For T = T c the negative eigenvalues of the Hessian become positive or zero. 
There is an exponentially large number of minima, connected by flat regions. 

• For T c > T > T K the number of relevant minima becomes proportional to 
exp(A r S*(/3)). where S*(/3) is the configurational entropy or complexity. The 
minima are separated by barriers that diverge with the number of particles 
N (in mean field theory), however the barriers are are finite in real life (i.e. 
beyond mean field theory). 

• For Tk > T the number of relevant minima is finite. The minima are sepa- 
rated by barriers that diverge with N. 

For T c > T > T K there is a dual description: the system may be described as 
a liquid and also as a solid (with an exponentially large number of different 
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solid structures). This duality tell us that all the information needed to study 
the glass can be extracted from the liquid phase. 



4 A more quantitative approach 

As we have seen we can write 

Z{(3) = £«P(-/W/«(/?)) = / dA^/,/3)exp(-/3A7), (6) 



where f a {(3) is the free energy density of the valley labeled by a at the tem- 
perature and Af(f,/3) is the number of valleys with free energy density 
less than /. 

We can assume that M(f ) f3) = exp(A r S(/, /?)), where the configurational en- 
tropy or complexity £(/,/?) is positive in the region f.fo(P) and vanishes 
at / = fo(/3). The quantity /o(/3) is the minimum value of the free energy: 
Af(f,P) is zero for f< f ((3) [6-9]. 



If the equation 



has a solution at / = f*((3) (obviously this may happens only for f*(f3) > 
/o (/?)), we stay in the liquid phase the free density is given by 

F = r-p- 1 nr,P) (») 



andE(r,/3) = E*(/3). 

Otherwise we stay in the glass phase and 

F = fM • (9) 

In order to compute the properties in the glass phase we need to know £(/, ($): 
a simple strategy is the following. 

We introduce the modified partition function 

Z( r ,P) =exp(-7V 7 $( 7;/ 3)) = £exp(- 7 N/ (/3)). (10) 
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Using standard thermodynamical arguments it can be easily proven that in 
the limit N — > oo one has: 



7«(T,«=7/-S0J,/), / = ^i«l. 



The complexity is obtained from $(7; /3) in the same way as the entropy is 
obtained from the usual free energy [12]: 

W/) = ^. (ID 



A few observations are in order: 

• In the new formalism 7, the free energy and the complexity play respec- 
tively the same role of f3, the internal energy and the entropy in the usual 
formalism. 

• In the new formalism (3 only indicates the value of the temperature which 
is used to compute the free energy and 7 controls which part of the free 
energy landscape is sampled. 

• When /3->oowe sample the energy landscape: 

Z(7;oo) = ^exp(-7iVe a ) = f u(e)deexp(-'jNe) (12) 

where e a are the minima of the Hamiltonian and v[e) the density of the 
minima of the Hamiltonian. 

• The equilibrium complexity is given by = S(/3; /3). On the other hand 
£(7; 00) give us information on the minima of the Hamiltonian. 

It is convenient to write 

Z( T , (3) = jdC exp (- 7 # (C) - N-ym C) + iV 7 /(7, G)) = 

J dCexp(-7//(C)-7V7/(/3,C) + 7V7/(7,C)), (13) 

where f(fl, C) = C) — /(oo, C), f((3, C) is is constant in each valley and 
it is equal to the free energy density of the valley to which the configuration C 
belongs, and C) is the temperature dependent part of free energy density 
of the valley to which the configuration C belongs. 

Our strategy is to use duality and liquid theory to evaluate the properties of 
f(P,C). 
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A simple approximation is the so called quenched approximation which in this 
case consists in assuming that 

eM-Af(P, C)) » exp(-A/ 7 (/3)), (14) 

where / 7 (/3) is the expectation value of f(/3,C) taken with the probability 
distribution proportional to exp(— ^yH(C)). 

The quenched approximation would be exact if all the valleys at a given tem- 
perature (also with different free energy) would have the same entropy. In other 
words the fluctuations of the energy are more relevant than the fluctuations 
of the entropy. 

It is now a matter of simple algebra to show that 

$(7; /3) = F liquid (i) + / 7 (/3) - / 7 ( 7 ), 
£(7; 0) = S liquid {l) - 5 7 ( 7 ) + / 7 (7) - W, (15) 

where = df 1 {(3)/d 1 . 

Everything is now computable using liquid theory. Different approximations 
can be done for / 7 (/3) and slightly different results may be obtained for the 
transition temperature and for the specific heat in the glass region. The replica 
formalism allow us to study directly the function $(7; (3) and most of the 
computation are done (for sake of compactness and simplicity) using this 
powerful tool. 



5 Some results 

Analytic and numerical computation have been done in many case. I will 
present here only the results for the Kob Andersen Lennard Jones binary 
mixture described in the introduction at density p = 1.2. 

In order to compute the free energy (or the entropy of the valleys) a simple 
and useful approximation is to assume that the potential near each minimum 
of the Hamiltonian is approximately an harmonic one and therefore the free 
energy of a valley can be computed using an harmonic approximation. In this 
way one finds that the previous formulae in the liquid phase become 

— Siiquid(P) — S harmonic {f3) + . . . , (16) 
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Fig. 4. Analytic and numerical complexity as function of T A . 

where 

^~ = ^ ndet (c4)^))' (17) 



Similar formulae holds below the transition temperature. 

The relevant quantities can be or computed analytically (using liquid theory) 
or extracted by numerical simulations in the liquid phase. 

Both the analytic and the numerical results of [9] for the configurational are 
shown in fig. 4 (from [9]). The correctness of the numerical estimate have been 
later confirmed by the more accurate computations of [13]. 

These microscopic analytic computations support the proposal that there is a 
thermodynamic liquid glass transition characterized by vanishing of the com- 
plexity and that at low temperature below Tr the system stays only in a few 
valleys. 

In other words in the glass phase replica symmetry is broken and we aspect 
that there are characteristic violations of the fluctuation dissipation theorem 
in off-equilibrium dynamics in the aging regime [14-16]. Indeed in the aging 
regime the fluctuation dissipation theorem cannot be applied and generalized 
dissipation relations are valid. The form of these generalized dissipation re- 
lations is fixed by the theory and the resulting predictions are confirmed by 
numerical simulations. 

Direct measurements of these fluctuation dissipation relations in real (not 
numerical [17]) experiments both for glasses are possible with the present 
technology. Some experiments are under the way. The outcome will be a crucial 
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test for this theoretical approach and I hope that some results will be available 
in the next future. 



References 



[1] For example see S.Sastry, P. Debenedetti, F.H.Stillinger, Nature 393 554 



(1998). K. K. Bhattacharya, K. Broderix, R. Kree, A. Zippelius, cond- 



mat/9903120| . L. Angelani, G. Parisi, G. Ruocco and G. Viliani, Phys. Rev. 



Lett. 81 4648 (1998) 

[2] A. Barrat, R. Burioni, and M. Mezard, J. Phys. A 29, L81 (1996), A. Barrat, 
S. Franz and G. Parisi J. Phys. A: Math. Gen. 30, 5593 (1997). 

[3] T.R. Kirkpatrick and D. Thirumalai, Transp. Theor. Stat. Phys. 24, 927 (1995) 
and references therein. 

[4] A. Crisanti and H.-J. Sommers, J. Phys. I (France) 5, 805 (1995). 

[5] M. Mezard, G. Parisi and M.A.Virasoro, Spin glass theory and beyond, World 
Scientific (Singapore 1987). 

[6] M. Mezard and G. Parisi, Phys. Rev. Lett. 82 747 (1998); A first principle 
computation of the thermodynamics of glasses x>nd-mat 9812180| , J. Chem. 
Phys. in press. 

[7] M. Mezard, How to compute the thermodynamics of a glass using a cloned 
liquid, |cond-mat/9812024| , J. Chem. Phys. m press. 

[8] B. Coluzzi, M. Mezard, G. Parisi and P. Verrocchio, Thermodynamics of binary 
mixture glasses, cond-mat/990312^ ., J. Chem. Phys. in press 

[9] B. Coluzzi, G. Parisi and P. Verrocchio, Lennard-Jones binary mixture: a 
thermodynamical approach to glass transition, |cond-mat / 9904124 , J. Chem. 
Phys, (in press); The thermodynamical liquid-glass transition in a Lennard- 
Jones binary mixture, |cond-mat /9906124 , Phys. Rev. Lett, (in press). 

[10] Silvio Franz, Giorgio Parisi Effective potential in glassy systems: theory and 
simulations, cond-mat/9711215. 

[11] Miguel Cardenas, Silvio Franz, Giorgio Parisi Constrained Boltzmann-Gibbs 
measures and effective potential for glasses . . . cond-mat / 9801 155 . 

[12] E. Monasson Phys. Rev. Lett. 75 (1995) 2847. 

[13] W. Kob, F. Sciortino, P. Tartaglia, Inherent Structure Entropy of Supercooled 



Liquids, x>nd-mat /9906081 



[14] L.F. Cugliandolo and J. Kurchan, Phys. Rev. Lett. 71, 173 (1993); J. Phys. 
A: Math. Gen. 27, 5749 (1994). 

[15] S. Franz and M. Mezard, Europhys. Lett. 26, 209 (1994) 



11 



[16] S. Franz, M. Mezard, G. Parisi and L. Peliti, Phys. Rev. Lett. 81 1758 (1998); 
The response of glassy systems to random perturbations: a bridge between 
equilibrium and off-equilibrium, |cond-mat /9903370| . 

[17] G. Parisi, Phys. Rev. Lett. 79, 3660 (1997). 



12 



